Scale dependence of halo and galaxy bias: effects in real space 
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We examine the scale dependence of dark matter halo and galaxy clustering on very large scales 
(0.01 < A;[feMpc~^] < 0.15), due to non-linear eflects from dynamics and halo bias. We pursue a two 
line offensive: high resolution numerical simulations are used to establish some old and some new 
results, and an analytic model is developed to understand their origins. Our simulations show: (i) 
that the 2 = dark matter power spectrum is suppressed relative to linear theory by 5% on scales 
0.05 < k[hMpc~^] < 0.075; (ii) that, indeed, halo bias is non-linear over the scales we probe and 
that the scale dependence is a strong function of halo mass. High mass haloes show no suppression 
of power on scales k < 0.07[/iMpc~^], and only show amplification on smaller scales, whereas low 
mass haloes show strong, ^ 5 — 10%, suppression over the range 0.05 < fc[/iMpc~^] < 0.15. These 
results were primarily established through the use of the cross-power spectrum of dark matter and 
haloes, which circumvents the thorny issue of shot-noise correction. The halo-halo power spectrum, 
however, is highly sensitive to the shot-noise correction; we show that halo exclusion effects make this 
sub-Poissonian and a new correction is presented. Our results have special relevance for studies of 
the baryon acoustic oscillation features in the halo power spectra. Non-linear mode-mode coupling: 
(i) damps these features on progressively larger scales as halo mass increases; (ii) produces small 
shifts in the positions of the peaks and troughs which depend on halo mass. We show that these 
effects on halo clustering are important over the redshift range relevant to such studies {0 < z < 2), 
and so will need to be accounted for when extracting information from precision measurements of 
galaxy clustering. 

Our analytic model is described in the language of the 'halo-model'. The halo-halo clustering term 
is propagated into the non-linear regime using '1-loop' perturbation theory and a non-linear halo bias 
model. Galaxies are then inserted into haloes through the Halo Occupation Distribution. We show 
that, with non-linear bias parameters derived from simulations, this model produces predictions that 
are qualitatively in agreement with our numerical results. We then use it to show that the power 
spectra of red and blue galaxies depend differently on scale, thus underscoring the fact that proper 
modeling of nonlinear bias parameters will be crucial to derive reliable cosmological constraints. In 
addition to showing that the bias on very large scales is not simply linear, the model also shows 
that the halo-halo and halo-dark matter spectra do not measure precisely the same thing. This 
complicates interpretation of clustering in terms of the stochasticity of bias. However, because the 
shot noise correction is non-trivial, evidence for this in the simulations is marginal. 

PACS numbers: 98.80.-k 



I. INTRODUCTION 



Statistical analysis of the large scale structures ob- 
served in galaxy surveys can provide a wealth of in- 
formation about the cosmological parameters, the un- 
derlying mass distribution and the initial conditions of 
the Universe. The information is commonly extracted 
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through measurement of the two-point correlation func- 
tion [ujsl or its Fourier space analogue the Power spec- 
trum J H) 0] ■ When further combined with high 
precision measurements of the temperature anisotropy 
spectrum from the Cosmic Microwave Background very 
strong constraints can be imposed on the initial condi- 
tions, the energy content, shape and evolution of the Uni- 
verse [§[. 

For homogeneous and isotropic Gaussian Random 
fields, such as is supposed for the post inflationary den- 
sity field of Cold Dark Matter (hereafter CDM) fiuctu- 
ations, each Fourier mode is independent, and thus all 
of the statistical properties of the field are governed by 
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the power spectrum. However, non-linear evolution of 
matter couples the Fourier modes together, and power 
is transfered from large to small scales 0, Conse- 
quently, it is non-trivial to relate the observed structures 
to the physics of the initial conditions. Further, since one 
typically measures not the mass, but the galaxy fluctu- 
ations, some understanding of the mapping from one to 
the other is required. This mapping, commonly referred 
to as galaxy bias, encodes the salient physics of galaxy 
formation. 

One last complication must be added: since galaxy po- 
sitions are inferred from recession velocities using Hub- 
ble's law, and because each galaxy possesses its own pe- 
culiar velocity relative to the expansion velocity, a non- 
trivial distortion is introduced to the clustering on all 
scales from the velocity field. These velocity effects are 
commonly referred to as redshift space distortions. Thus 
one must accurately account for non-linear evolution of 
matter fluctuations, bias and redshift space distortions 
in order to extract precise information from large scale 
structure surveys and this remains one of the grand chal- 
lenges for modern physical cosmology. 

In this paper, we investigate the issue of bias in some 
detail, through both numerical and analytic means. We 
focus on real space effects and reserve our results from 
redshift space for a subsequent paper. Our numerical 
work focuses on the generation of multiple realizations of 
the same cosmological model, in two different box sizes. 
This allows us to construct halo catalogues spanning a 
large dynamic range in mass that are largely free from 
discreteness fluctuations and the multiple realizations al- 
low us to derive errors that are 'true errors from the en- 
semble'. We use this data to show that not only is halo 
clustering on very large scales scale dependent, but that 
the scale dependence is a strong function of halo mass. 
These results are completely expected given the standard 
theoretical understanding of dark matter haloes based on 
the 'peak-background split' argument [Til, fl^. fisl. [l^. [Tsl . 

Our analytic approach to modelling these trends can 
be summarized as follows. 

• Haloes are biased tracers of the mass distribution. 
To describe this bias, we assume that the bias re- 
lation between the halo density field and the dark 
matter field is non-linear, local and deterministic. 
This allows us to use the formalism of [H, [l^ . In 
order for the local model to hold, one must inte- 
grate out small scales where locality is almost cer- 
tainly violated. This can be done in real space by 
smoothing small-scale fluctutations, or in Fourier 
space by considering small wavenumbers |^ . 

• The underlying CDM density field is then prop- 
agated into the non-linear regime using standard 
Eulerian perturbation theory techniques [lol| . 

• Galaxies are then assumed to form only in haloes 
above a given mass [s^ and these are inserted into 
each halo using the 'halo model' approach [l^, [26l . 



Because we write the PT evolved halo density field as 
a series expansion we refer to this method as 'Halo-PT' 
theory. 

Our results are particularly relevant for studies which 
intend to use the baryon acoustic oscillation feature 
(hereafter BAO) in the low redshift clustering of galaxies 
to derive constraints on the dark energy equation of state 
0. The CDM transfer function that we have adopted 
throughout contains a significant amount of BAOs, and 
we give an accounting of the possible non-linear correc- 
tions from mass evolution and biasing that might influ- 
ence the detection and interpretation of such features. 
Previous work in this direction has primarily focused on 
analysis of numerical simulations [Sg, [13, [H, [s^, H^l , al- 
though several analytic works have recently been pre- 
sented: [11] derive the exact damping of BAOs in the 
Zel'dovich approximation and calculate it in the exact 
dynamics by re-summing perturbation theory; [4^ con- 
sider real space corrections to the power spectrum from 
one-loop PT; [ij use the halo model, also in real space, 
to explore systematics; [i^ [1^ consider a model of La- 
grangian displacements fit to simulations. 

In Section [IT] we discuss how the approach we have de- 
veloped here is complementary to and expands on these 
studies. In Section lllll we discuss the numerical simula- 
tions and present our measurements of scale dependence 
in the dark matter, halo center and halo-dark matter 
cross power spectra. We also present the evidence for 
large scale non-linear bias. Then in Section llVi we out- 
line some key notions concerning the halo model of large 
scale structure, as this is the frame work within which 
we work. In Section [IV CI we describe the non-linear bias 
model that we employ. In Section llVDI we use the 3rd Or- 
der Eulerian perturbation theory to describe the evolved 
Eulerian density field in terms of the inital Lagrangian 
fluctuations. In Section|Vl we use the 3rd order halo den- 
sity fields to produce an anlytic model for the '1-Loop' 
halo and halo-cross dark matter power spectra. In Sec- 
tion [VTj we explore the predictions of the analytic model 
for a range of different halo masses. In Section IVIII we 
compare our analytic model to the nonlinear bias seen in 
the numerical simulations. We use the analytic model to 
examine the galaxy power spectrum in Section IVini and 
present our conclusions in Section 1X1 

Throughout, we assume a flat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) cosmological model with en- 
ergy density at late times dominated by a cosmological 
constant (A) and a sea of collisionless cold dark mat- 
ter particles as the dominant mass density. We take 
rim = 0.27 and Ha = 0.73, where these are the ratios of 
the energy density in matter and a cosmological constant 
to the critical density, respectively. We use a linear the- 
ory power spectrum generated from cmbfast [47], with 
baryon content of fli, ~ 0.046 and h = 0.72. The nor- 
malization of fluctuations is set through trg — 0.9, which 
is the initial value of the r.m.s. variance of fluctuations 
in spheres of comoving radius 8h~^ Mpc extrapolated to 
z = using linear theory. 
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II. MOTIVATION 

A number of recent papers have attempted to quan- 
tify the scale dependence of galaxy bias. A subset of 
these have forwarded simple analytic models to remove 
the scale-dependent biases in the power spectrum esti- 
mator. We discuss some of these below so as to set the 
stage for our work. 



A. Cole et al. (2005) 

Based on analysis of mock galaxy catalogues from the 
Hubble volume simulation, these authors proposed a sim- 
ple analytic model to account for the non-linear scale de- 
pendence: 



1 



1 + Aik 



(1) 



The parameter Ai = 1.4, and A2 was allowed to vary 
over a narrow range, which was then marginalized-over 
in the fitting procedure. When Aik ^ 1, this model has 



A.k^] 



(2) 



We show below that the bracketed terms are suggestive 
of the Pii+Pii terms from the dark matter perturbation 
theory, but with incorrect dependence on fc. In addition, 
ignoring the fact that Ai may depend on galaxy type is 
a serious inconsistency. For instance, our results indicate 
that Al for LRG-like galaxies is smaller than 1.4. 

A further concern regarding this model is that no ac- 
counting for non-Poisson shot noise has been made. If 
galaxy formation takes place only in haloes, then galax- 
ies are not Poison samples of the mass distribution. The 
analytic model we develop shows that it is important to 
account for this, and how. We are therefore skeptical 
about the blind use of equation ([1]), particularly with 
regard to its use in the analysis of LRGs [sO] . 



B. Seo &c Eisenstein (2005) 

These authors examined the scale dependence of halo 
bias in a large ensemble of low-resolution numerical sim- 
ulations. They proposed 



(3) 



which bares some similarity to that of Cole et al., as 
can be seen by re-writing Ai — > Ai/(6^PLin) and A2 — > 
A2 / [h^ P-Lin) ■ The effective spectral index of the linear 
power spectrum evolves from > rteff > —1 over the 
range of k of interest, so we can think of Plih as being 
approximately constant. Then, the inclusion of the 6^ 
term decreases the effective Ai as required, but it also 
decreases the effective A2] our results indicate that this 
is inappropriate. 



There is an important difference between this model 
and the previous one — the inclusion of the constant 
power term Aq. This was introduced to account for 
'anomalous power' - by which was meant effects envis- 
aged by [29j - and or non-Poisson shot noise following 
|48|. Our analyis strongly supports the inclusion of this 
term. 



C. Seljak (2001); Schulz & White (2006); 
Guzik, Bernstein & Smith (2006) 

These authors explored the scale dependence of galaxy 
bias in the halo model focusing their attention on the 1- 
Halo term. They showed that if this was taken simply as 
a non-Poisson shot-noise correction, then it would be a 
significant source of scale dependence in the large scale 
galaxy power spectrum. We have confirmed this in our 
study. These studies lead one to suggest a base form of 
the kind: 



Pgg(fc) =62pLin(fc)+Ao 



(4) 



On top of this base form we need to include modifications 
that capture the true non-linear evolution of the halo 
field. 



D. Huff et al. (2006) 

These authors used a set of three large cosmological 
simulations to investigate the scale dependence of the 
BAOs. They suggested 



Pgg(fc) = fe'FLi„(fc)exp -(Aifc) 



An 



(5) 



where the exponential damping term was introduced to 
account for halo profiles. Again examining the large scale 
behaviour, Aik <^ 1, we see that this equation may be 
re- written 



Pss{k)^b''Pun[l~{Aikf] +Ao. 



(6) 



If < ^1 < 1 (the range they considered), then this 
formula will suppress the power spectrum on scales k < 
0.1/iMpc^^ by a percent at most. If this model is to 
account for the non-linear corrections that we see, then 
Al > 1. However, the strong exponential damping makes 
it unlikely that this model will properly characterize the 
transition from the 2- to 1-Halo term. This is because, 
as we show below, the nonlinear evolution of halo centers 
includes an additional boost at intermediate k which this 
model does not capture. 



E. A necessary model 

Our results suggest that a necessary model will have 
the following properties: the model should be able to pro- 
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duce a pre-virialization feature and a small scale nonlin- 
ear boost with fc-dependencies motivated by physical ar- 
guments; nonlinear corrections should depend on galaxy 
type; a constant power term should be added to account 
for non-Poisson shot noise. We therefore expect a rea- 
sonable starting point for any empirical modelling of the 
large scale, scale dependence of the galaxy power spec- 
trum to be, 



x\W{k/ko.)\ + Ao{T)\W{k/kp)\ + — .(7) 

rig 

Our notation makes explicit that the coefficients have the 
following properties: b > 0, Aq > 0, Ai > 0, A2 > and 
m > 0, and that all depend on galaxy type T. The first 
term is composed of two pieces: in the first piece, we have 
modeled the damping of BAOs using a Gaussian as de- 
rived in [41I, \^ for the dark matter case; for the second 
piece, we have added a fc-dependent boost that models 
the power added by mode-mode coupling and nonlinear 
bias. Our simple power-law form (with two parameters 
A2 , m) is meant to describe this effect over a restricted 
range of scales. The fact that this term is additive as op- 
posed to multiplicative, is meant to emulate the fact that 
in PT this corresponds to the term which arises from 
the convolution of linear power on different scales and 
is therefore smooth possesing no information on BAOs. 
For weakly nonlinear scales it has a positive spectral in- 
dex (note, however, that in the limit fc — > 0, m = 4 is 
expected from momentum conservation arguments) . The 
second term corresponds to Poisson shot noise from un- 
equal weighting of haloes. The last term corresponds to 
the Poisson shot noise from the galaxy point distribution. 
We have included filter terms W{k/ka) and W{k/kp) to 
indicate the damping due to density profiles, which will 
occur for k > l/rvir- This function may be greatly sim- 
plified by examinig the case k ^ l[/iMpc~^], for which it 
reduces to 

Pgg(fc,T) = b\T){Pun[l-Ai{T)e]+A2{T)k"'^^^} 

+Ao{T), (fc« 1) , (8) 

where the parameter Aq (T) subsumes all sources of con- 
stant large scale power. 



III. SCALE DEPENDENT HALO BIAS FROM 
NUMERICAL SIMULATIONS 

A. Simulation details and halo catalogues 

We have performed a series of high-resolution, coUi- 
sionless dark matter iV-body simulations, where N — 
512^ equal mass particles. Each simulation was per- 
formed using Gadget2 51]; the internal parameter set- 
tings can be found in Table 1 of [52], where more de- 
tails about the runs themselves are available. The initial 



TABLE L Halo samples. A'^icai is the number of independent 
realizations. iVh and nh are the ensemble average number and 
number densities of haloes in each mass bin. 









L[Mpch~^] 


iVh [xlO''] 


fih [Mpc~^/i^] 


LR 


Bin 1 " 


8 


1024 


3.6863 


3.43313 X 10"'' 


LR 


Bin 2 ' 


8 


1024 


7.3530 


6.8480 X 10"^ 


HR 


Bin 3 " 


20 


512 


6.9287 


5.1623 X 10""* 


HR 


Bin 4 


20 


512 


5.5415 


4.1287 X 10"^ 



"Mass bin 1 = Af > 1.0 X lO^^h-^ ATq 

Mass bin 2 = 1.0 X lO^^h-'^MQ > M > 4.0 X IO^^/i-^Mq 
'^Mass bin 3 = 4.0 X lO'^'-^h''^ Mq > M > 7.0 X lO^'^h'^MQ 

Mass bin 4 = 7.0 X IOI^/j-iMq > M > 4.0 x lO^^/i-iM, 
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conditions were set-up using the 2nd-order Lagrangian 
Perturbation theory at redshift z = 49 [s^, with linear 
theory power spectrum taken from cmbf ast }47j , with the 
cosmological model being the same as that used through- 
out this paper. We will present results from two differ- 
ent box sizes: 20 smaller higher resolution box (hereafter 
HR) for which the volume is V = = (512/i-iMpc)3, 
and 8 realizations of a larger, lower resolution box (here- 
after LR) for which L = 1024ft." ^Mpc box. 

Haloes were identified in the z — Q outputs using the 
friends-of-friends algorithm with linking-length parame- 
ter I — 0.2. Halo masses were corrected for the error 
introduced by discretization of the halo density struc- 
ture [5^. Since the error in the estimate of the halo 
mass diverges as the number of particles sampling the 
density field decreases, we only study haloes containing 
50 particles or more. For the HR and LR simulations 
this corresponds to haloes with M > 4.0 x lO-^^ft.~^M0 
and M > 4.0 x IO^^H'^Mq, respectively. We then con- 
structed four non-overlapping sub-samples of haloes with 
roughly equal numbers per sub-sample. The two low 
mass bins were harvested from the HR simulations and 
the two high mass ones were taken from the LR runs. 
Further details may be found in Table HI 



B. Mass, halo and halo-Mass power spectra 

For each realization and each bin in halo mass we mea- 
sured the following quantities: the power spectrum of the 
dark matter P^^{k); the power spectrum of dark mat- 
ter haloes P^^{k); and the cross power spectrum of dark 
matter and dark matter haloes, P^^{k). The power spec- 
tra may be generally defined 

((5"(k),5'3(k')) = {2TT)H^{k + k')P"^(fc) , (9) 



where 



P'\k) 



p^^{k) 



P^^ik) 
P^'^ik) 



(10) 

with 5{k) and ^'^(k) being the Fourier transforms of the 
mass and halo density perturbation fields. 



p'(x) = p-[l + <5«(x)] 



(11) 
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0.05 0.1 0.15 0.05 0.1 0.15 



k [h/Mpc] k [h/Mpc] 




0.05 0.1 0.15 0.05 0.1 0.15 

k [h/Mpc] k [h/Mpc] 



FIG. 1: Halo power spectrum measurements and predictions, ratioed with a smooth 'No Baryon' dark matter power spectrum, 
for four bins in halo mass: results for massive haloes (the top two panels) are from the LR simulations, whereas lower masses 
(the two bottom panels) are from the HR simulations. Filled circles in the top, middle and bottom sections of each panel show 
the ensemble average non- linear P^\k) P'^'^(fe), and P^^, respectively. The open circles in the middle sections show P^^{k) 
with the non-standard shot noise subtraction described in Appendix |A] In all panels the linear theory dark matter, halo-halo 
and halo-cross power spectra are shown as dashed lines. The top panel also shows predictions from halof it (solid lines) [s^ ] 
and 1-Loop perturbation theory (dot-dash lines). Solid lines in the middle and bottom panels show our new analytic model, 
Halo-PT. 



6 



where the index i again distinguishes between dark mat- lowing estimator for the halo bias matrix 



ter and haloes, e.g. = p and p^ = p^. 

We estimate the spectra through the conventional Fast 
Fourier Transform (FFT) method [95^ (for a detailed dis- 
cussion see [54j). The mean power and 1-a errors on the 
spectra were estimated from the ensemble neglecting the 
bin-to-bin covariances. Inspection of an estimate of the 
covariance matrix from the 20 HR simulations showed 
that this is reasonable. There is, however, a small degree 
of off diagonal covariance, but the number of simulations 
was insufficient to make a precise estimate. 

Figure [1] shows the three types of power spectra mea- 
sured from the ensemble of simulations for each of the 
four bins in halo mass described in Table ID The filled 
circles and associated error bars in the top, middle and 
bottom sections of each panel show P^^{k), P^'^(fc|Af), 
and P^^{k\M). The open circles show the result of apply- 
ing a non-standard shot- noise correction to P^^, which 
we describe in Appendix [Al 

To emphasize the non-linear evolution of the spectra 
and the BAOs in particular, we have divided each spec- 
trum by a smooth linear theory spectrum, which we shall 
refer to as our 'No Baryon' model. This was constructed 
by performing a chi-squared fit of the cmbf ast transfer 
function data to the smooth transfer function model of 



Tik)=\l + {aq+ibq)'-'^ + {cqfy 



-l/d 



(12) 



The derived parameters are: q = fc/0.19, a — 4.86, b — 
4.81, c= 1.72, d= 1.18. 

To compare the halo spectra with the linear theory 
we require estimates of the halo bias on very large scales, 
which we measure as follows. We begin by assuming that 
the spectra can be written 



P'^{k\M)^A'^Pun{k) ; A'' 



1 b^^ 



(13) 



where i and j denote the type of spectrum considered 
and where, for reasons that will later be apparent, we 
distinguish between the bias from P^^ and P^^. Hence, 
the likelihood of obtaining an estimate of the power P*-' 
in the Z*^ fc-bin is assumed to be an independent Gaussian 
with dispersion ai 



C{P;',ai\A,Pun) = 



1 



exp ■ 



2af 



(14) 

Thus the combined likelihood of obtaining the data set 
, (T;| can be written. 



CaPn,{<yi}\A,Pun) = n ^Pr,<^l\APLin) . (15) 



1 = 1 



A^ 



1 = 1 



Puniki)Pl' /af 



(16) 



We construct error estimates through further differenti- 
ation of the Gaussian likelihood function: 



9^ log C 
~dAP~ 



-1/2 



-1/2 



On maximizing the likelihood function, we find the fol- 



(17) 

Lastly, since we observe scale-dependent non-linear ef- 
fects in the matter power spectrum for k > 0.05/iMpc~^, 
we only use modes with k < 0.04/i Mpc^^ in the fitting 
of the amplitude matrix A. Estimates of the large scale 
bias parameters b^^ and b^^ are presented in Table [TTl 

M 

The dark matter power spectra in Fig. [T] show sig- 
nificant deviations away from linear theory prediction: 
at 0.05 < k[hMpc^^] < 0.075, there is a suppres- 
sion of power relative to linear theory, whereas at fc > 
0.075 [/iMpc~^] there is an amplification. Perturbation 
theory studies [l^l refer to the suppression effect (caused 
by tidal terms) as 'pre-virialization'. Recently, this has 
been understood in much more detail as a result of the 
damping of linear features by nonlinearities, leading to 
an exponentially decaying propagator that measures the 
loss of memory of the density field to the initial condi- 
tions [111, \^ . Although the effect in the power is rather 
well-known and has been observed in recent numerical 
simulations of CDM spectral^, HI, HI, H^l, our results 
constitute a rather precise measurement of this effect, 
with realistic errors drawn from the ensemble. A com- 
plete assessment of the damping of linear features such as 
BAOs is done by studying the propagator [4ll . [43 | rather 
than the power spectrum. See (ssj for further discussion 
on this. 

The solid and dot-dashed lines show predictions based 
on halof it 54] and on 1-Loop perturbation theory (PT). 
The PT results do well compared to the simulations on 
very large scales, but for k > 0.07 they increasingly over- 
predict the power. We note also that PT appears to 
predict the pre-virialization feature in the simulations, 
adding additional support to the claim that P^^{k) re- 
quires non- linear corrections on the scales of interest. 
However, qualitatively, it under-predicts the magnitude 
of the effect. We note that halof it does reasonably well 
at capturing the behavior for k < 0.07hMpc~^, but it 
appears to under predict the measured data. 

Before moving on to P^^^ and we think it is worth 
noting that the BAOs in P^^ on fc > 0.1/iMpc"^ have 
been erased. The large-box LR measurements show that 
the third peak is gone, and the height of the second 
peak has dropped so that it appears more as a plateau. 
However, the behavior at fc < 0.05ft.Mpc~'^ appears un- 
changed. 

Discreteness corrections for P^^{k) have been studied 
in some depth[5J|. However, for P^^, the appropriate 
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correction is more complicated because haloes are rare, 
highly clustered, and spatially exclusive. In Appendix 
|A]we show that the standard Poisson shot noise correc- 
tion for the cluster power spectrum results in negative 
power at high k. This lead us to propose a new method 
for making the 'shot-noise' correction that accounts for 
exclusion, which we discuss therein. The open circles 
in the middle sections of each panel in Figure [T] show 
the result of this new correction. Filled circles show the 
uncorrected power, and stars show the standard correc- 
tion — clearly, the choice of correction is crucial. Note 
that owing to the arbitrary normalization things for the 
standard shot noise method look better than they actu- 
ally are. Unfortunately, the residual uncertainties in our 
new procedure prevent us from making strong statements 
about the scale dependence of halo-halo clustering. 

Whilst the discreteness correction is troublesome for 
P^^ it is almost negligible for P^^ (the halo model ar- 
guments which follow allow us to quantify this). Our 
estimates of P^^ are shown in the bottom sections of the 
panels in Fig[TJ Notice that the scale dependence of P^^, 
is a strong function of halo mass. P^^ for the most mas- 
sive haloes shows no deviations from linear theory until 
k > 0.7/iMpc~^. However, the pre-virialization feature 
appears and gets enhanced as one goes to lower masses. 
Indeed, for our lowest mass bin, P^^ is sub-linear until 
k > 0.15/iMpc~^ 

This has important consequences for the BAOs. In the 
highest mass bin (top left panel), the oscillations in P^^ 
at fc > 0.07/iMpc~^ have been erased. However, the first 
trough, at A: ~ 0.04/iMpc~^, is unaffected. For the next 
mass bin (top right panel) , the first peak and trough are 
unmodified, and that the second peak is becoming notice- 
able. This trend continues as we decrease mass; there 
is even a hint of the third peak in the bottom panels. 
These measurements indicate that non-linear dynamics 
can erase oscillations on progressively larger scales as halo 
mass increases and that small displacements to the posi- 
tions of the peaks and troughs may occur and that these 
will also be dependent on halo mass. If the locations of 
these peaks and troughs are to play an important role 
in constraining cosmological parameters, our measure- 
ments suggest that understanding and quantifying these 
displacements will be very important. 

Before continuing, we comment on the possible expla- 
nation of these shifts through simple scatter from cosmic 
variance. We remark that it is certainly possible to rec- 
oncile some of the shifts in the peak positions through 
this. However, we draw attention to the fact that all of 
the points k > 0.05 in the cross-power spectra of the low 
mass haloes are systematically lower than expected from 
the linear theory. We also re-iterate that the derived er- 
ror bars are the errors on the means for 20 realizations. 
One caveat is that since the spectra are normalized by 
the very-large scale modes of the power spectrum, where 
cosmic variance errors are larger, we expect some small 
fluctuations in the relative amplitudes of the theory pre- 
dictions as more data is acquired. Estimates of the error 



TABLE II: Bias parameters for the halo samples, bi , and 
bf^ are the first three non-linear halo bias parameters derived 
from the Sheth & Tormen model p^ . [l5l | averaged over halo 
bins, bi, ti2 and 63 are the parameters measured from the 
scatter plots, b^^ is the large scale bias parameter measured 
directly from P*'\ 

b'^ tr br br tj bj &r 

Bin 1 " 2.28 ±0.03 2.19 0.94 -0.98 2.23 1.68 -4.08 

Bin 2 1.49 ±0.04 1.53 -0.30 0.94 1.41 0.04 -1.29 

Bin 3 1.02 ±0.03 1.13 -0.46 1.47 1.04 -0.85 0.37 

Bin 4 0.87 ± 0.03 0.98 -0.44 1.44 0.91 -0.74 0.55 

"See Table |I] for definition of bins. 

in the present LR simulations suggest changes of the or- 
der ^ 1% to be acceptable; and this increases to ~ 3% 
for the HR simulations. If the amplitudes for the theory 
curves are too low by ~ 3%, then some of these discrep- 
ancies may be alleviated. However, it is unquestionable 
that non-linear effects are present on these scales and we 
must therefore firmly accept that it is likely that these 
may cause some shifting of the harmonic series. Only 
a wider and expanded numerical study will be able to 
address and answer these questions more completely. 

The solid lines in the middle and bottom sections of 
each panel show predictions from the analytic model de- 
scribed in the following sections. In all cases this model 
provides a better description than does linear theory. 

C. Scale dependence of the bias 

Next, we examine the scale dependence of halo bias. 
We will consider 

?;hh /L ^ / P^^(fci) t5\iii,\ P^^jkj) 

ONL(fc») = \ - , ; °Nl(«*) = ^7777- • (18) 

as well as 

For any particular realization the wave modes of the halo 
and dark matter density fields arc almost perfectly corre- 
lated. Because the first set of estimators are derived from 
taking the ratio of measured power spectra, they are in- 
sensitive to this source of cosmic variance. In this sense, 
the second set of estimators are non-optimal. However, 
they are the ones which will be used with real data, since 
P^^ is generally not observable. 

In Figure [5] we show the results of measuring these 
quantities for the same halo mass bins as in the previous 
section. The top and bottom parts of each panel show 
(fT5)) and (HH]), respectively. The error bars, which were 
derived from the ensemble to ensemble variations, are sig- 
nificantly larger for than for (|18p. as expected. The 
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FIG. 2: Large scale halo bias derived directly from A'^-body simulations for four bins in halo mass. The top sections of each 
panel show the estimators bf^ (solid blue points) and 6nl with and without shot noise correction (large and small red stars, 
respectively). The bottom panels show the same, but for bf^^ and bj^^ (See equations [18] and [19] for definitions). The solid blue 
and red lines in each panel show the predictions for the bias from our Halo-PT model. 
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solid lines show the predictions from the new analytic 
model described in the next section. 

For b^^ we show both the shot-noise corrected (large 
stars) and un-corrected (small stars) results. As was the 
case for the halo-halo power spectra, we see that this 
correction is important, so it must be known rather ac- 
curately. The estimators for b^^ are also shown (filled 
circles). Except for the highest mass bin, b^^ > b^^. 
Indeed, as we shall argue later, there are compelling the- 
oretical reasons why the biases derived from P^^ and P^^ 
are not in fact the same, and that, one generally expects 
^hh > ^Sh^ However owing to the uncertainty regarding 
the shot noise correction, no firm statement can be drawn 
from the current data. 

Two possible explanations why the highest mass halo 
bin appears to behave differently are: Firstly, if the shot 
noise correction to P^^ is too aggressive, we may have un- 
derestimated the true P^^ and therefore the bias. Alter- 
natively, we have made no shot-noise correction to P^^; 
if one should be applied (we think this is unlikely, owing 
to the large number of dark matter particles), then our 
current estimate of b^^ may be biased high. 

The estimators based on equation (fT9|) show more scale 
dependence than those based on (fTHl) . especially for the 
two high mass bins. As noted by [4j, this is because the 
BAOs are erased on larger and larger scales as higher and 
higher mass haloes are considered. Thus on dividing the 
halo spectra by a linear theory BAO spectrum, we are in 
fact introducing scale dependence from the linear model. 

For the highest mass haloes, b^]^ is constant at fc < 
0.04/iMpc"\ but it increases monotonically as fc in- 
creases. This is a direct consequence of the absence of 
pre-virialization in P^^ on intermediate scales and the 
rapid onset of non- linear power on smaller scales, com- 
pared to P^^ . The bias bfl\^ is flat for haloes with masses 
(4.0 X 10^3 < M[h-^MQ] < 1.0 X 10^'^), suggesting that 
small clusters and group mass haloes are linearly biased 
tracers of the non-linear dark matter. The smallest two 
bins in halo mass show the reverse trend: the bias de- 
creases at large fc, by ~ 10% compared to the approxi- 
mately constant value at smaller fc. 



from some mass function and with the density pro- 
file of each halo being drawn from some universal 
stochastic profile. 

The model attains its full potential when the second as- 
sumption is stated: 

• All galaxies exist only in isolated dark matter 
haloes, with more massive haloes hosting multiple 
galaxies. 

In essence the model has been in existence for several 
decades [H, M, M, EH, [Hi, [el, [H. However it was 
not until the advent of large numerical TV-body simula- 
tions and accurate characterization of halo phenomenol- 
ogy that its true value was realized. Namely, given an 
appropriate Halo Occupation Distribution (HOD - i.e., 
a prescription for the number and spatial distribution of 
galaxies within a halo) the model successfully reproduces 
the real-space form of the two-point correlation function 
of galaxies over a wide range of scales. It predicts subtle 
deviations from a power-law which have recently been 
seen in observations [1^ and provides a framework for 
describing the luminosity ^od, and environmental de- 
pendence of galaxy clustering 68, 69] . It also enables new 
tests of the CDM paradigm to be constructed [tO, [Zlj . 



B. Power spectrum 

In the model the density fields of haloes and dark mat- 
ter may be written as a sums over haloes, 



E 



(20) 



where i ~ {1,2} distinguishes between haloes and dark 
matter and Ni = N and N2 = iVh are the total number of 
haloes and the number of haloes in some restricted range 
in mass. M, and Xj are the mass and center of mass of 
the jth halo and Uj = pj{xj)/Mj is the mass normalized 
density profile. Following Ji^J, the power spectra P^^ , 
P^^ and P^^ can be written as the sum of two terms: 



IV. THEORETICAL MODEL 



P'^ik) 



(21) 



We now describe a model for interpreting the trends 
seen in the previous section. Our discussion is based 
on the 'halo-model', which we briefly summarize below. 
See [591 for a more detailed review. 



The first term, P^^' referred to as the '1-Halo' term, 
describes the intra-clustering of dark matter particles 
within single haloes; the second, -P2H' referred to as the 
'2-Halo' term, describes the clustering of particles in dis- 
tinct haloes. They have the explicit forms: 



A. The 'Halo Model' of large scale structure 

The halo model may be described by the simple state- 
ment: 

• All dark matter in the Universe is contained within 
a distribution of CDM haloes, with masses drawn 



^Mk) = 



1 



dMn{M)M^ |C/(k|M)|^ 



x0y(M, M) 
2 

Yl {dMin{Mi)MiUi{k\Mi)} 



(22) 



1 



P'P^ Jo 



1=1 



X P^'' {k\Mi , M2 ) e,, {Ml , M2 ) , (23) 
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where n{M) is the halo mass function, which gives the 
number density of haloes with masses in the range M to 
M + dM, per unit mass. The 6ij matrix carves out the 
halo density field to be considered, e.g. for haloes with 
mass M > Afcut the matrix is 

e,,(Mi,M2) = 

e(Mi)e(M2) e(Mi)e(M2 - Meut) 

_ e{Mi - Meut)e(M2) e(Afi - Meut)e(M2 - Meut) _ 

where <d{x) is the Heaviside Step Function. More com- 
plicated halo selections can easily be described through 
the 9y notation. Lastly, P^^{}<i\Mi, M2) is the power 
spectrum of halo centers with masses Mi and M2 ■ This 
function contains all of the information for the inter- 
clustering of haloes; precise knowledge of this term is 
required to make accurate predictions on large scales. 

In principle, P^^{k.\Mi,M2) is a complicated func- 
tion of Ml, M2 and k. Initial formulations of the halo 
model (isl I27I [H, It^I assumed that it could be well- 
approximated by 

Pi^\fc|Mi, M2) = 6i(Mi)&i(A^2) PLin(fc) , (24) 

where all the scale dependence is in PLin(fc), which is 
taken from linear theory, and all the mass dependence is 
in the scale-independent bias parameter hi (M) [l^, 13 ■ 
In this approximation, Pe^'^(k|Mi, A/2) is a separable 
function of Mi, M2 and k. As we show below, compar- 
ison with numerical simulations shows that this simple 
model over-predicts power on very large scales and pro- 
vides insufficient power on intermediate scales. In both 
cases, this is about a ten percent effect. 

This discrepancy is not unexpected [H, [tI, [tII- A 
simple correction results from setting 

phh(fc|Mi, M2) = bi{Mi)hi{M2)PMk) , (25) 

where Pnl is the non-linear rather than the linear matter 
power spectrum, and, in addition, imposing an exclusion 
constraint: 

4"'(r|Mi, M2) = -1 ; (r < r^-^ + ^vir J , (26) 

where rvir is the virial radius of the halo. 

The success of this approach is demonstrated by com- 
paring P^^{k) measured in the z = output of the Hub- 
ble Volume simulation [t^ with the halo model calcu- 
lation. The open and filled symbols in the top panel 
of Figure [3] show the measurement before and after sub- 
tracting a Poisson shot noise term (which is shown by the 
triple dot-dashed line.) The dot-short dash line shows the 
linear theory prediction, and the other two dot-dashed 
curves show two estimates of the 2-halo term: the one 
which drops more sharply at large k is based on equa- 
tion ((25)) and the other one is based on the original ap- 
proximation of equation (j24p . Since equation (|25p re- 
quires the use of a nonlinear power spectrum, we used 
the one provided by [53 |. 



The symbols in the bottom panel show the measure- 
ments divided by the halo model calculation which uses 
equation ([M)) for the 2-Halo term (i.e. the initial linear 
theory-based approximation). Notice how they drop be- 
low unity at fc ^ 0.1/iMpc~^. The solid line shows the 
halo model calculation which is based on equation — 
it reproduces this pre-virialization feature well. 

As an interesting aside, we note that the fitting formula 
halof it does very well at matching the pre-virialization 
feature. Whilst it is not apparent from the figure we 
also point out that the transfer function of the Hubble 
volume simulation does contain BAOs; thus, our results 
demonstrate that the fitting formula of [3 appears to 
be accurate, for this data, for BAO models to roughly 
~ 5%. In hght of this, we note that the discrepancy 
between halof it and the mass power spectra from our 
smaller box HR simulations is somewhat puzzling. We 
highlight this issue for further study, one possible expla- 
nation is the difference in initial power spectra, on the 
other hand, also note that a calculation within the frame- 
work of Renormalized Perturbation Theory [4l| suggests 
that even for these simulation volumes on expects small 
effects due to absence of coupling to large scales. 

To fix the small discrepancies which remain, some 
authors have advocated making the halo bias factors 
scale dependent [t^, but the implementation has been 
based on fitting formulae rather than fundamental the- 
ory. While equation (P5)) appears to fare better than the 
original approximation (j24p , as we will soon show, in go- 
ing from to ([^5]) , one is making the assumption that 
halo bias is linear even when the mass density field is not. 
If this is not the case then the method is incorrect. In 
addition, there is an unpleasant circularity in requiring 
prior knowledge of P^i^{k) in order to predict PNL(fc). 



C. Halo bias: The non-linear local bias model 

The discussion above makes clear that a rigorous treat- 
ment of the 2-Halo term is currently lacking. This term 
requires a description of how dark matter haloes cluster. 
Whereas current models seek to describe halo clustering 
as a biased version of dark matter clustering, the scale 
dependence of halo bias is still rather poorly understood 

m 112, [il Q H Izi H H H m M, Hi- in the 

following sections we develop a model to understand its 
main properties. In particular, we will discuss a gen- 
eral non-linear, deterministic, local bias model for dark 
matter haloes. This model is exactly analogous to that 
derived for galaxy biasing by [T^ and first applied to dark 
matter haloes by [l3| . 

To begin, consider the density field of all haloes with 
masses in the range M to M -I- dM, smoothed with some 
filter of scale R. We now assume that this field can be 
related to the underlying dark matter field, smoothed 
with the same filter, through some deterministic mapping 
and that this mapping should apply independently of the 
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FIG. 3: Dark matter power spectrum measured as a function 
of wavenumber measured from the 2 = time slice of the 
Hubble Volume simulation [t^. In the top panel the red 
and blue points show the estimates of the dark matter power 
spectrum measured from the simulation, with and without a 
Poisson shot noise correction. The dot-dash line shows the 
linear theory and the triple dot-dash line shows the Poisson 
correction. The blue dotted and dot-dash curves show the 
1- and 2-IIalo terms. The thick red dot-dash curve shows 
the 2-IIalo term where Pnl has been used instead of Pliii. 
The bottom panel presents the ratio with respect to the halo 
model, but with equation (|24p for the 2-Halo term. The solid 
red line shows the effect of the Pnl modification. 



precise position, x, in the field: i.e. 



(27) 



where the sub-scripts on the function J-' indicate that it 
depends on the mass of the haloes considered and the 
chosen filter scale. The filtered density field is 

^(x|i?) = d\S{y)Wi\^-y\,R) , (28) 

VF(|x|,i?) being some normalized filter. Taylor expand- 
ing TiAf^ii} about the point S ~ yields 



h{M\R) 



[S{^\R)Y 



(29) 



We now assume that there is a certain filter scale above 
which J^{M,R} is independent of both the scale considered 
and also the exact shape of the filter function. Hence, 



S^\x\R,M) 



1=0 



where the bias coefficients are 



(30) 



(31) 



The linear bias model has 6^ = for alH > 1. 

The bias coefficients from the Taylor series are not in- 
dependent, but obey two constraints. The first arises 
from the fact that ((5^(x|i?)) = 0, which leads to 



(32) 



Thus, in general 6o is non-vanishing and depends on the 
hierarchy of moments. This allows us to re- write equation 

as 

i=l 

(33) 

Nevertheless, we may remove 6o from further consider- 
ation by transforming to the Fourier domain, where it 
only contributes to (5(k = 0). 

The second constraint states that a sum over all halo 
density fields 5^{x\M) weighted by halo mass and abun- 
dance must recover the dark matter density field p^ . 
This requires that 



1 



dMn{M)Mbi{M) 



1 {^ 
{i 



1) 

0,2,3,, 



■) 



(34) 



For CDM models whose initial density perturbations 
are Gaussian Random, the bias coefficients may either 
be derived directly through the 'peak-background split' 
argument [ll|, [H, [H, [l3| or measured directly from N- 
body simulations. Figure [4] shows the halo bias param- 
eters up to third order, derived in the context of the 
Sheth-Tormen mass function; see [iH] for the analytic 
expressions. We compare these with measurements from 
our simulations in Appendix IB] 

A practical application of this method rests squarely 
upon our ability to truncate the Taylor series at some 
particular order. However, since the procedure that we 
have adopted for doing this requires some further knowl- 
edge, we shall reserve our discussion until Section IVB I 



D. HaloPT: Evolution of halo fields 

We now evolve the halo density field as expressed by 
equation (|30p into the non-linear regime via perturbation 
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FIG. 4; First three halo bias parameters derived from the 
Sheth-Tormen !14| mass function as a function of halo mass 
[l5| . The solid line shows bi, the dashed line 62, and the 
dot-dashed line 63. 



theory techniques. For a short discussion of these meth- 
ods see Appendix [Cl and for a full and detailed review 
see [lo| . The main idea that we require from perturba- 
tion theory is that each Fourier mode of the density field 
may be expanded as a series, 



^(k,a)=^i?l(a)(5,(k) 



(35) 



where (5i(k) is the ith order Eulerian perturbation and 
Di{a) is the linear growth factor. Thus, on Fourier trans- 
forming the halo bias relation of equation truncated 
at third order, and on inserting the FT expansion from 
above, we arrive at (keeping up to cubic terms) 



<5\k|M,i?) - 6i(A/) [(5i(k|i?) + <52(k|i?)+J3(k|i?) 
^ 62 (Af) f d\ 



2! 



^i(qi|i?)(5i(k-qi|i?) 



(27r)3 

+ 2,5i(q|i?)^2(k-qi|i?) 

J {2nf ^Mi\R)5i{ci2\R) 
X 5i(k-qi-q2|i?) , (36) 



<5h(k,a|M,i?) = ^i?5'(a) [<5h(k|M)], 



(37) 



n=l 



(27r)3""3 
xF„h(qi,...,q„|M,i?) 



(38) 



where [(5''(k|Af )]„ is the nth order perturbation to the 
halo density field, and where the short-hand notation 

[5^(k)]^^ = (5'°(k - qi q„) has been used. The 

functions i^J^ (qi, q„|A/, i?) are the Halo-PT kernels, 
symmetrized in all of their arguments. The first three 
may be written in terms of the dark matter FT kernels: 



= h^{M)W{\k\R)F^ ■ 

FI2 = h^{M)W{\k\R)F^,2 
, h2{M) 



W^(|qi|i?)W^(|q2|i?)FiF2 



(39) 



(40) 



^^^2,3 = h^{M)W{\k\R)F^a,3 



62 (M) 



W^(|qi+q2|i?) 



-l¥(|q3|i?)Fi,2F3-f2cyc 
63 (M) 
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W{\cii\R)W{\ci2\R)W{\ci:i\R)FiF2F3 (41) 



where F^ ,^ = i^j' (qi, q^ |M, i?). Thus equations 
(157)1 and ([38l) can be used to describe the mildly non- 
linear evolution of dark matter halo density fields to ar- 
bitrary order in the dark matter perturbation, and equa- 
tions (P5VHT|) make explicit the halo evolution up to 3rd 
order. Together, these ideas define our meaning of the 
term 'Halo-FT'. 

It is now apparent that halo clustering studies which 
assume a linear bias model and take the power spectrum 
to be the fully non-linear one, are effectively assuming 
that 5''(k|Af) = &i(M)X;'5j(k), with hi{M) < 61 (M) 
for i 7^ 1. However, for CDM, the peak-background split 
argumeirt informs us that this never happens, unless the 
density field itself is linear [111, IHj [H) Hj, [3 • We must 
therefore conclude that extrapolating the linear bias rela- 
tion into the weakly non-linear regime, without full con- 
sideration of the non-linearity of halo bias is incorrect. 



V. THE 1-LOOP HALO MODEL 



A. Halo center power spectra 



where 5i(q|i?) = W{\c{\R)5i{q). We next insert the solu- 
tions for each order of perturbation, which are presented 
in equation (jC4|) of the Appendix, into equation ([36)) . 
On re-arranging terms and collecting powers of 5, the 
mildly non-linear density field of dark matter haloes may 
be written as a FT series expansion of the dark matter 



We now use the Halo-FT to calculate the power spec- 
trum of halo centers in the mildly non-linear regime. We 
define the power spectrum of halo centers for haloes with 
masses Mi and M2 to be 

(J,\ki|Afi)'5c(k2|A/2)) = (27r)35^(ki2)Pc'^ki|Mi,M2) . 

(42) 
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On inserting the Halo-PT solutions for each order of the 
perturbation we find that P^^'^{k) can be written as the 
sum of three terms 

P^''{k\A'h,M,) = Pl%{k\M,,M2) + Pl%{k\M,,M2) 
+2P^%{k\M,,Nh) , (43) 

where 

phh^(fc|Mi,M2) = F^{k\Nh)F^{k\M2)Pii{k) ■ (44) 

P,%(A;|Mi,Af2) = 2 j ^Pn(9)Pii(|k-q|) 

X F2^q, k - q|Afi)P2h(q, k - qjMa) ; (45) 

P,%(fc|Mi,M2) = 6Pn(fc) j ■j0^Piiiq) 

X F^^{k,q\M,,M2) . (46) 

Here Pii(fc) is equivalent to the linear theory power spec- 
trum and 

P^f(k,q|Afi,M2) - i [P^(k|Mi)P3^(k,qi,-qi|M2) 
+pi^(k|M2)P3^(k,qi,-qi|A/i)] . (47) 
I 



When these expressions are averaged over all halo masses, 
weighted by the respective cosmic abundances Min{Mi) 
and M2n{M2), then the constraint equation p4|) guaran- 
tees that they reduce to the standard '1-loo p' e xpression 
for the PT power spectrum of dark matter [1C) |: 



PlLoop (fc) = Pll (fc) + P22 (fc) + Pi3 (fc) . (48) 



Strictly speaking the 1-Loop power spectrum refers to 
P22 + P13, we shall break convention and use equation 
(US)) to define what we mean. Explicit details of the 1- 
Loop expressions may be found in Appendix IC 21 



The theory may be further developed by directly sub- 
stituting the Halo-PT kernels, given by equations ([35]- 
BT|) . into equations piHi^ . A little algebra shows that 



P^''{k\M,,M2,R) 



6l(A/i)6i(Af2)PlLoop(fc) + 6l,3(Afi, A//2)c7^(P)Pii(fc) 



\WikR)\" 

+2 6i,2(Mi,Af2) J ^l^l^^M|pMp^^(^){p^^(|k-q|)P2(q,k-q) + 2Pn(fc)P2(k,-q)} 
h2{Mi)b2{M2) f d\ |W^(gP)|^ |T4^(|k-q|P)| 



d\ W^(gP)VK(|k-q|P) 

Pii(q)Pii(|k-q|) , (49) 



(2^)3 \W{kR)\ 



2 



where 

6,.,(Afi, Af2) = \ [b,{Mi)bj{M2) + b,{MiMM2)] . 

(50) 

Before continuing, we point out and answer an important 
question that naturally arises at this junction: How does 
one compare the filtered theory with the unfiltered obser- 
vations? We forward the proposition that the unfiltered 
nonlinear power spectrum can be recovered through the 
following simple operation: 



This is unquestionably true for an observed non-linear 
field. It will therefore also be true for the correct theo- 
retical model. 

This rather lengthy expression may be more readily 
digested through the examination of two limiting cases. 
But first, notice the important fact that it is still a sepa- 
rable product of mass dependent terms and scale depen- 
dent terms. 



P{k\R)/\W{kR)\^ ^ P{k) . (51) 



When the two halo masses are identical, then 
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+bi{M)b2{M) 



Wm? Ai(9)Ai(|k-q|) . (52) 



This expression is equivalent to evolving the non-linear, 
local, galaxy bias model [lB| through Eulerian PT. This 
has been explored by [111 and [sst - 

Secondly, consider the case where we integrate over 
one of the halo masses, say M2, weighting by M2 and its 



abundance n{M2). Equation again insures that all 
terms involving b2{M2) and 63(^1^2) vanish, and so the 
resulting expression is the 1-Loop correction to the halo 
center-dark matter cross power spectrum: 



\W{kR)\^ 



fci(Af)PiLoop(fc) + -53(A/)a"(i?)Pn(fc) \ +b2{M) 



(Pq W{qR)WlSk-ci\R) 



W{kR) 



X Pn(q) Fn(|k-q|)F2(q,k-q) + 2Pn(fc)^^2(k,-q) 



(53) 



Inspection of these two limiting cases reveals three re- 
markable features: 

• Firstly, if the non- linear bias parameters 62 (M) 
and 63 (M) are non- vanishing then the bias on large 
scales is not &i(M). 

• Secondly, the halo-halo spectrum has a term that 
corresponds to constant power on very large scales, 
whereas the cross spectrum does not. 

Both of these points were independently noted by 
I 



and [S^l, but for the case of non- linear galaxy biasing 
(also see nil). 

• Thirdly, the large scale bias derived from the halo- 
dark matter cross power spectrum is not 61, nor 
is it given by the bias derived from the halo-halo 
power spectrum. 

To see these points more clearly we take the k ^ 
limit of equations (gH]), JS^l) and ((55|) : 



Pr(fc|Afi,M2,i?) = 6i(Mi)6i(Af2)Pii(fc) |1 + a'{R) 

, 62(Mi)52(M2) f d\ 



34 (h2{Mi) , b2{M2 



(27r)= 



21 \hi{Ah) bi{M2) 
Pii{q)\W{qRt'^ 



1 A3(Mi) 63(^2) ^ 

2 ^6i(Mi) bi{M2), 



(54) 



P^'^mM.R) = biiM)P,iik) \ l + a\R) 



68 b2{M) bsjM) 
21 61 (M) ^ bi{M)_ 



bl{M) 



^[Pi,{q)\W{qR)\']' ■ (55) 



P°'^{k\M,R) = bi{M)Pu{k) { 1 + a^{R) 



M b2{M) I bsjM) 
216i(Af) 2 6i(Af) 



(56) 



These expressions make the first point noted above trivially obvious: The large scale bias is modulated by 
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the Halo-PT correction terms, and these depend on the 
non-Unear bias parameters 62 and 63 and also on the fil- 
tered variance of fluctuations. 

The second point noted above originates specifically 
from the quadratic non-linear bias terms found in equa- 
tions and ([SS]) . e.g. terms containing b2{Mi)b2{M2) 
and b\{M). For a linear power spectrum that obeys the 
limit Pii(fc ^ 0) — > 0, these expressions reduce to the 
constant 



P^^{k 0|M) 



hl{M) f d\ 



[Pn{q)\WiqR)\' 



(57) 

This term was discussed in great detail for the case of 
galaxy biasing by [23 |. 

The third point noted above can be understood by con- 
structing the Hnear bias, e.g. dividing equations ([56]) and 
(155)1 by Pii. On squaring the bias recovered from ([5^ 
and subtracting it from the bias from (j55p . we find 



2Pn{k)J (2^)3 



[Pn{q)\W{qR)\' 



-a\R) 



^62(M) + i63(M) 



(58) 
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FIG. 5: '1-Loop' bias parameters l}^"^ and h^"^ in the ultra-large 
scale limit. The solid through to dotted curves show 6*''^, mea- 
sured on scales (feobs = {0.001, 0.005, 0.01, 0.05} /iMpc-^), 
derived from the 1-Loop halo-halo cross power spectrum as 
a function of halo mass. The triple-dot dash curve shows 
6**^ as derived from the 1-Loop halo-dark matter cross power 
spectrum. 



B. Convergence of the power spectrum 



We now see that, because P approaches a constant 
on very large scales, on dividing through by Pii(fc|A/) 
the bias function 6'^'^ diverges at the origin as 1/P^Y^ 
diverges. 

Figure [5] shows our expressions for and feLhi (from 
equations [54] and [55]) . In this particular case, we as- 
sume the nonlinear bias parameters derived from the 
Sheth-Tormen mass function by [l^. We use a Gaus- 
sian filter, for which W{kRQ) = exp[— (fci?)^/2]. For 
Rg = 20/i-iMpc and for our fiducial cosmology, we find 
a'^iRc = 20) « 0.046. To inspect the differences more 
closely, we take the ratio of the predictions with respect 
to the tree-level theory, e.g. hi = P{^/Pii ^/P^l'/Pn. 
The figure demonstrates two of the points raised above. 
Firstly the large scale bias is not simply bi : halo-halo bias 
(solid through to dotted curves) does not converge as one 
considers larger and larger scales; however bias from the 
cross power spectrum is very close to linear for all ex- 
cept the most massive haloes, where 62 and 63 are very 
strongly rising functions (see Fig. . Secondly, it is now 
obvious that 6'^'^ and b^^^ are not the same. Note also that 
the magnitude of the expected scale dependence: the bias 
varies by at most five percent when k is changed by an 
order of magnitude. The mass dependence of b^^ shown 
in Fig. [5] is simply driven by that of 63 (M) (again see 
Fig. |4]) . As halo mass increases the bias slowly increases 
until it reaches a maximum at M '--^ 10^^h~^MQ. It then 
decreases to M ~ 10^'^h~^MQ after which it shoots up 
dramatically for larger masses. 



We now return to the issue of truncation and applica- 
bility of the Taylor series expansion of the halo field. A 
first requirement for the Taylor series to converge after 
a finite number of terms is that the filter scale be large 
enough so that the rms dark matter fiuctuations be much 
less than unity: (t{R) <C 1. Consider the case where R 
is very large and convergence occurs at first order, we 
then have: S^{x\R) = bi{M)Si{x\R). As the filter scale 
R is slowly decreased the rms fluctuations in 6 increase 
and a larger and larger number of terms are required to 
accurately map the underlying bias function. Finally, as 
a{R) 1 all terms in the series are required. At this 
point the method has no merit. 

Since a robust criterion for truncation is out of reach at 
present, we propose an ad-hoc criterion for convergence 
that must plausibly be obeyed, that is 

bi{M)a{R) < 1 . (59) 

In Appendix [B1 we shall also discuss an empirical method 
for testing convergence. 



C. Returning to the halo model 

We now translate these ideas back into the language of 
the halo model. To begin, we shall restrict our attention 
to the halo model in the large scale limit, more precisely 
we consider scales where U{k\M) ~ 1. Since the Fourier 
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transform of the mass normalized profile may be written 



U{k\M) = / d\U{r\M) 



sin[A;r 
kr 



(60) 



our large scale condition simply becomes fcrvir ^ 1- 
In the above equation we have, for convenience, as- 
sumed spherical density profiles. If halo mass and virial 
radius are related through M = 47r200pr^ij./3, then 
for the largest collapsed objects in the Universe M « 
10^^/i~^Mo, and the above condition translates to the 
inequality k <^ 0.4/iMpc~^. If we assume an NFW den- 
sity profile [11] with concentration parameter c ~ 6, then 
for k ~ 0.15/iMpc"\ we find that U{k) « 0.994. We 
therefore assume that this is an excellent approximation 
over the scales that we are interested in. 

Hence for scales fcrvir ^ 1 , our equations and 
at the 1-Loop level in Halo-PT, now take the forms: 



\W{kR)f 



|W^(fci?)|^ 



1 



P'P' Jo 
xQ,,{M,M) ; 

1 



dMn{M)M'^ 



(61) 



P"P^ Jo 



Y[{dMiniMi)Mi} 



1=1 



phh(k|Afi,M2,fl) 

\W{kR)f 



ey (Mi,M2) , (62) 



where we have explicitly included a filter on the 1-Halo 
term. (Recall that the halo center power spectrum in the 
2- Halo term already includes such a filter.) 

We now see that, because the bi{M) are the only mass 
dependent functions, on insertion of equation (|49|) into 
equation (j62p the integrals over mass may be immediately 
computed. Thus, 



(63) 



where P'^^'^^{'k\R) is equivalent to equation (|49p. except 
that we have replaced all of the mass dependent bias 
parameters by the average ones: i.e. 

^ dMn{M)Mb,{M)Q{M - Afeut) 
^ J^' dMn{M)Me{M - Mcut) ' 

We note that if we wish to weight by halo number density 
rather than mass density then we simply remove the mass 
weighting in the numerator and denominator of (bi). 



VI. EVALUATION OF THE THEORY 

In this section we present the results from the direct 
computation of the 1-Loop halo center expressions for 
P^^ and P^^, as given by equations ([5^ and re- 
spectively. This section is almost entirely pedagogical; 
we urge those who are only interested in the direct com- 
parison with the numerical work to press on to Section 

Em 



Recall that it is necessary to adopt some filter scale R. 
We have studied two choices: Rq = 20ft.~^Mpc and Rq = 
lO/i^^Mpc for which the linear theory, Gaussian filtered, 
variances are <t^{R) — 0.046 and 0.177, respectively. The 
larger smoothing scale is required for the more massive 
haloes. 



A. Halo— dark matter cross power spectra 

In Figure [6] we show the predictions for the scale de- 
pendence of P*^ at the 1-Loop level, as a function of wave 
number. The different panels show results for different 
halo masses and smoothing scales. In all four panels, 
we see that, as expected, there is a small (few percent) 
positive offset from the linear theory bias value bi. The 
largest offset occurs for the cluster mass haloes, but here 
it may be the case that, owing to the bias being large, the 
filter scale that we have adopted for these objects may 
still be too small for adequate convergence. We also note 
that the offset for the M = lO^^ft-~^M0 haloes is negligi- 
ble. This can be attributed to the fact that 62 ~ &3 ~ 
(see Fig. |4]). For the lower mass haloes the offsets are 
roughly ~ 2% in excess of linear. 

Considering the predictions for the highest mass 
haloes, we see that the spectrum is scale independent 
up to A: ~ 0.07/iMpc~^, where the non-linear amplifi- 
cation from the P22 term becomes dominant. For this 
case, the absence of the pre-virialization feature may be 
understood as follows: Firstly, we note that 61 and 62 are 
positive, whereas 63 is negative. However, since 63 is an 
order of magnitude smaller than the others, it plays no 
significant part in determining the shape of the spectrum. 
The P13 term in the 1-Loop power spectrum, which is the 
main cause of the pre-virialization feature, is thus over- 
whelmed by the action of the quadratic bias 62 term. 
This suggests that P^^ should be scale independent up 
to fc ~ 0.07/iMpc"^ 

Next, we collectively consider the predictions for the 
lower mass haloes, as they show many similar traits. 
Firstly, we find that when k < 0.01/iMpc~^, then the ra- 
tios of the 1-Loop to tree level (or linear) spectra are flat. 
However, for O.Ol/iMpc"^ < k < 0.07/iMpc"\ significant 
scale dependence is apparent: the pre-virialization fea- 
ture is present and it appears to become stronger as halo 
mass decreases. On smaller scales still, the non-linear 
boost from the P22 term amplifies the power spectrum 
and breaks all scale independence. Interestingly, the on- 
set of P22 is pushed to smaller scales as halo mass de- 
creases. These effects can be understood as follows. For 
these objects bi and 63 are positive, whereas 62 is neg- 
ative. On large scales, we see that 62 and 63 are nearly 
equivalent, but 63 is slightly dominant, and this results 
in a small positive correction. Whereas on smaller scales 
this trend reverses and 62 becomes dominant. The overall 
correction is then negative and this leads to the enhanced 
pre-virialization feature and delay of the onset of P22 . 

It is also interesting to note the imprint of the BAO fea- 
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FIG. 6: 



The real space cross power spectrum of haloes and cold dark matter at the 1-Loop level on large scales. 

.15 1 nl4 1 nl3 



The four 



panels show predictions for haloes with masses: M — {lO^^, 10^'', 10^^, 10^^} Mq. In the upper plot of each panel, the 
thick (black) curve represents the total spectrum as given by equation (|53p . The dash curve gives the 1-Loop contribution from 
the linear bias parameter bi term; the dotted curve and dot-dash curves give the 1-Loop contributions from the non-linear bias 
parameters 62 and 63, and the thickness (color) of the lines indicates their sign, with thick (red) lines being positive and thin 
(blue) lines being negative contributions. The thin solid (green) line gives the smoothed, linearly biased, linear power spectrum. 
The lower plot of each panel shows the ratio of the total cross spectrum with the smoothed linear theory cross spectrum. 
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FIG. 7: The real space halo-halo power spectrum at the 1-Loop level on large scales. The four panels show predictions for haloes 
with masses: M — {10^^, 10^*, 10^"^, 10^^} h~^MQ. In the upper panels of each sub-figure, the thick (black) curve represents 
the total spectrum as given by equation (I52|l . the dash curve represents the 1-Loop contribution from the pure linear bias (6i) 
term, the dotted and dot-dash curves represent the 1-Loop contributions from the &261 and 6361 terms, respectively. The triple 
dot-dash curve represents the quadratic nonlinear bias term h\. Lastly, the thin solid (green) line denotes the smoothed, linearly 
biased linear theory spectrum. As before, line thickness (color) distinguishes positive and negative contributions. Bottom panels 
show the ratio with the linear spectrum. 
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tures in the ratios of the power spectra. The strength of 
the signal appears to depend on halo mass and increases 
as halo mass decreases. As we discussed in Section ITlI C| 
this can be attributed to the fact that non-linear evolu- 
tion suppresses BAOs on small scales, and on taking the 
ratio with a linear theory spectrum, we are artificially 
introducing oscillations. 

B. Halo-Halo power spectra 

In Figure [7| we show the predictions for the scale de- 
pendence of P^^, at the 1-Loop level, as a function of 
wavenumber. The four panels show results for haloes 
with masses in the same range as Figure [51 Again, upper 
panels show the contributions from each of the Halo-PT 
terms in equation (j52p . and where the sign of each con- 
tribution is distinguished through line thickness/color. 
Sub-panels are as before. 

Some obvious similarities exist between the auto-halo 
spectra and the cross spectra. In particular: the large 
scale bias is not given by bi ; the highest mass halo spec- 
trum shows no sign of the pre-virialization power decre- 
ment; the non-linear boost occurs increasingly at smaller 
scales as halo mass decreases; with the exception of the 
highest mass haloes, the ratio of the 1-Loop spectra to 
the tree level spectra have BAOs imprinted. These effects 
may all be understood through the explanations from the 
previous sub-section. 

We also notice some important differences between P^^^ 
and P^^. Firstly, the addition of the quadratic non-linear 
bias term, [62 (M)]^, modifies the results on the largest 
scales. As was discussed in Section |Vl ^'I'^ioop becomes 
a white-noise power spectrum as Pn 0, unless 62 — 
0. In the figure, the contributions from this term are 
denoted by the triple dot dash hues. On considering all 
four panels and paying special attention to the ratios, we 
see that, with the exception of the case M ^ 10^'^h~^MQ, 
there is an upturn in power on scales k < O.Ol/iMpc""'^. 
For the M 10"/i"^Mq haloes this effect is not found, 
this owes to the fact that 62 ~ (see Fig S]) . 

Secondly, we note that on smaller scales the effect 
of the 62 term is to boost the power across all scales. 
Thus the pre-virialization power decrement is no longer a 
decrement relative to the linear theory. However, relative 
to the power measured on say k ~ 0.01/iMpc~^, there 
is a very real decrement. Moreover, the fc-dependence of 
the terms multiplying 6| will make the decrement appear 
larger than we would expect from the non-linear matter 
spectrum. 



C. Non-linear evolution of BAOs 

We now briefly consider how mode-mode coupling and 
non-linear biasing affect the evolution of the BAOs. Fig- 
ure [5] compares the 1-Loop predictions for P^^ with the 
linear theory predictions. We again consider the set 



of halo masses: M = {lO^^, lO", 10", lO^^} /i-iMg, 
and, to emphasize the evolution of the BAO features, 
we have taken the ratio of each spectrum with the 'No 
Baryon' model of equation (|12p . In addition, because 
we are purely interested in the scale dependence, we re- 
normalize all ratios to unity at fc ~ O.Ol/iMpc"^. 

The top left panel shows results for the highest mass 
haloes: all but the first of the BAO troughs has been re- 
moved through non-linear evolution. This is a much more 
aggressive non-linear evolution than one expects from 
considerations of the dark matter alone. Secondly, the 
trough appears to have been displaced towards lower fre- 
quencies. These effects arise because both the quadratic 
bias (6|) and 6261 terms are positive, and so dominate 
over the negative 63 and P^^ terms. This fact, coupled 
with the injection of power from P22, means that the 
non- linear boost occurs at fc ~ 0.05/iMpc^^. Thus, the 
overall effect on the halo spectrum is to shift the first 
trough to smaller k. 

All the spectra of the lower mass haloes display a 
pre-virialization feature, the strength of which increases 
as halo mass decreases. In addition, the number of 
peaks and troughs which remain in the evolved spec- 
tra increases as halo mass decreases, because the strong 
non-linear amplification from the P22 term is delayed by 
the negative 62^1 corrections. However, in contrast to 
the high mass haloes, for the low mass haloes the first 
acoustic trough is shifted towards larger k because, when 
0.01 < k < 0.05, then the negative correction from the 
6261 term is dominant and so subtracts power from the 
higher frequency side of the oscillation. This acts to shift 
the overall pattern to higher frequencies. 

VII. COMPARISON WITH SIMULATIONS 

In this Section, we compare the predictions from the 
theoretical model with the results from the numerical 
simulations presented in Section IIIII 



A. Halo center power spectra re-visited 

Returning to our analysis of Figure [U We are now in 
a position to comment on the success of the Halo-PT 
model in comparison to the numerical simulations. The 
analytic model was evaluated using the semi-empirical 
bias parameters that were determined as described in 
Appendix |B] Since here we are concerned purely with 
the scale-dependence of the spectra and not their overall 
amplitude, we allowed the very large scale normalizations 
to be considered as free parameters. These were fit-for in 
exactly the same way as was done for the linear theory 
models in equation . Whilst this re-scaling would not 
be necessary if the bias parameters that we were adopt- 
ing were precise and accurate, since we have not been 
able to establish this in a robust manner we feel that this 
approach is acceptable, given that it has been equally ap- 
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FIG. 8: Non-Linear evolution of the BAOs as traced by the halo-halo power spectrum for four different halo masses. The solid 
(red) line denotes the the non- linear (1-Loop) spectrum as given by equation (|52|) : the dash (blue) line denotes the linear theory 
(tree level) power spectrum. We have taken the ratio of each spectrum with a smooth 'No Baryon' model as described in the 
text, and have renormalized each spectrum so that they are unity on a scale k ~ 0.01/iMpc~^. 



plied to the linear theory. Moreover, this is the method 
of analyzing real data. 

These predictions are shown as the solid lines in Fig- 
ure [TJ Comparison with the measured P^^, shows that 
the model and the shot-noise corrected data show rea- 
sonable correspondence over scales k < 0.07/i~^Mpc. On 
smaller scales the match is poor. However, the uncer- 
tainty in the shot-noise correction makes it very difficult 
to draw firm conclusions. 

Considering P^^'^, we find that the analytic model is in 
good agreement with the simulations over a much larger 
range of scales. The model correctly captures the mass 
dependence of the pre-virialization feature ~ low mass 
haloes show greater loss of memory to the initial den- 
sity fluctuations. The shifting of the nonlinear boost 
to smaller scales as halo mass decreases is also matched 
rather well. 

To assess whether -Pc^^iLoop ^ better fit to the sim- 
ulation data than is linear theory, we have performed a 
likelihood ratio test, assuming that the likelihood func- 
tions are Gaussian. In this case, a necessary statistic for 
model selection is for 

LK = rj-^ > 1 . (65) 

where the subscript 'max' refers to the parameter choices 
in the models that maximize the likelihood. Re- 
stricting the information to be A: < 0.1/iMpc^^, we 



find, going from low to high mass haloes, that LR = 
{1.06, 1.03, 1.05, 1.07}, respectively. 



B. Halo bias re- visited 

We now examine how well the Halo-PT model does 
at matching the nonlinear scale dependent bias of the 
halo centers as measured in the numerical simulations 
(Figure [2]) . The top section of each panel shows that the 
analytic model (solid blue and red lines correspond to 6^*1, 
and b^, respectively) captures, qualitatively, the scale 
dependence of the bias. The model shows a bias that 
increases with k for the high mass haloes but decreases 
with k for lower masses. However, there are some notable 
discrepancies: The model under predicts and then over 
shoots the measured relationship for the most massive 
haloes; for the next bin in halo mass (top-right panel), 
the measured bias is flat, whereas the model predicts a 
down turn after k ~ 0.1/iMpc^^; the model fares better 
in the two lowest mass bins, but the down-turn at high 
k is not seen in the data. However, we must stress that, 
with the exception of the M ^ 10^'^h~^MQ haloes, the 
model out-performs linear theory, which would predict 
constant bias on all scales. 

Having extolled the virtues of our model we now draw 
attention to its short-comings. Whilst the predictions 
provide a good match to the halo power spectra, they 
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do not simultaneously provide a good match to the scale 
dependence of the bias. If we re-consider our measure- 
ments of the non-linear matter power spectrum (upper 
sections in Fig [1]) P^^ , we see that the 1-Loop model 
(dot-dash lines) over predicts the LR simulations on 
scales k > 0.05ft.Mpc~^ and the HR simulations on scales 
k > O.OT/iMpc""'^. It is therefore unlikely that the model 
as presented here can be made to work precisely . 



VIII. GALAXY POWER SPECTRUM 



m Section Eel it is a very good approximation to set 
U^{k\M) = 1 when kr^i^ < 1. Thirdly, the constant 
pre-f actors transform as l/p^ — > where 



j dMn{M) {N\M) e(M - 



And finally, the galaxy bias parameters are 
1 



(67) 



dMn{M){N(M))h(M) . (68) 



How does the scale dependence of the bias depend on 
galaxy type? The answer has important consequences 
for future galaxy surveys that will measure the cluster- 
ing of specific sub-classes of objects. We may address 
this question using the halo model by changing the mass 
weighting in the integrals to: 



n{M)M 
n{M)M^ 



n{M) {N\M) 
n{M) {N{N - 



1)|M) 



(66) 



where {N\M) and {N{N - l)\M) are the first two facto- 
rial moments of the halo occupation probability function 
P{N\M), which gives the probability for a halo of mass 
M to host N galaxies. Use of the factorial moments of 
P{N\M) subtracts off a term corresponding to the self- 
correlation of galaxies e.g. f(r ^ 0) = 6^{r)/n. This 
corresponds to the Poisson shot noise term in Fourier 
space. Secondly, the mean density profile of dark mat- 
ter is changed to the mean density profile of galaxies in 
the halo: U{k\M) Us{k\M). Following the discussion 



These changes in equations (|6T|) and ([62]) yield the 1- 
Loop halo-model prediction for the galaxy power spec- 
trum. Explicitly: 



\W{kR)\^ 



\W(kR)\^ 



(69) 
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where again we have explicitly included a filter function 
on the 1-Halo term. On inserting our expression for P^^ 
from equation (|52p into the 2-IIalo term, and again notic- 
ing that in the large scale limit the mass integrals may 
be performed directly, we find: 



2H 



\W{kR)\' 



{ [bf]'PiLoop{k)+blbla'{R)P,iik)} 



^Pn(.)Pn(|k-q|)»ffll^i:i^ 



(27r) 



\WikR)\' 



W{qR)W{\k-q\R) 
WikR) 



{Pn(|k-q|)f^2(q,k-q) + 2Fn(fc)F2(k,-q)} , (71) 



r 



The resulting expression for the galaxy power spectrum 
is identical to that given by [13, [s^ , with one important 
difference — we include the large scale constant power 
originating from Pf^{k\R). 

To study the expected differences between red and blue 
galaxies, we use the parametric forms for (N(M)) mea- 
sured by ^] in the semi-analytic models of [Sq : 



(iVslAf) =0.7 



M 
Ms 



{Nr\M) 



(72) 

The blue galaxy parameters are: Mb = 4 x 10^^ h ^Mq; 
for haloes with masses in the range 10^^ < M/h~^MQ < 
4.0 X 10^^ then as = 0.0, for larger mass haloes as = 0.8. 
The red galaxy parameters are: = 2.5 x lO^^/i^^Af©; 



for haloes with masses greater than the cut-off mass an = 
0.9. For the second moment of the HOD we follow the 
model of Kravtsov et al. [Hi, so that {N{N - 1)|M) = 
{N\Mf - 1. This makes P{N\M) sub-Poissonian as sug- 
gested from the observations J23 [67I [7^ and the semi- 
analytic models [1^, [2^, m, l88l . 89j |: and secondly, this 
choice allows the first moment alone to fully specify the 
hierarchy of moments. 

In practice, we use the models above but impose 
a lower mass cutoff of Mmin = lO^^/i^^Af©. This 
yields ub = 4.10 x lO'^/i^Mpc"^, ur = 7.93 x 
10-3/i3Mpc-3, {bf = 1.20, 6f = -0.14, 6f = 1.06} and 
{6f = 1.39, = 0.09, b§ = 0.89}. (This estimate of 6f 
agrees with the observational determinations of simi- 
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FIG. 9: Scale dependence of the power spectrum of blue and red galaxies (thick solid blue and red lines) on very large scales. 
Top left panel shows the results obtained from the evaluation of equations (|69|l and (|71[1 . The 1-Halo term is represented by the 
thin solid lines. For the 2-Halo term we have the following contributions: terms denoted by dash lines; 6261 terms denoted 
by dotted lines; 63&1 terms denoted by dot-dash lines; 62 terms denoted by triple-dot-dash lines. Note again for the 2-Halo 
contributions, line thickness dictates sign. The middle panel shows the ratio of the red and blue galaxy 2-Halo terms (thick 
red and thin blue lines, respectively) with the no-Baryon linear model of equation (|12|l . The dash line shows the linear theory 
model with BAOs. The bottom panel shows the same as the middle only this time we show the effect of including the 1-Halo 
term. The right hand panel is similar to the left, only we exchange PiLoop in the 2-Halo term for the halofit model [sj. 



lar galaxies from the PSCz survey [33].) These values 
are easily understood by noting how bi(M) depends on 
halo mass (e.g. Figure 2]), the weightings given in equa- 
tion ([7^ . and recalling that the halo mass function de- 
clines exponentially at M > 10^^h~^MQ. Notice that 
the red galaxy bias parameters are all positive, whereas 
6^ for the blue galaxies is negative. Therefore, while we 
expect to see a pre-virialization feature in Pb(/c), we do 
not for the red galaxies. 

The left panel of Figure [S] shows the power spectrum 
of blue and red galaxies evaluated using the 1-Loop Halo 
model. The top section shows the individual contribu- 
tions from all the linear and nonlinear terms. Note that, 
for the blue galaxies the non-linear correction terms 6f 6f 
(blue dotted lines) and bfb^ (blue dot-dash Hues) are 
roughly the same order of magnitude as the 1-Halo term 
(thin blue line). However, for the red galaxies, the 1-Halo 
term (thin red line) dominates over the non-linear bias 



corrections by factors of a few. For both populations the 
quadratic bias terms [6|]^ (triple-dot-dash lines) appear 
to be negligible. 

The middle section of the left hand panel shows the ra- 
tios of the 1-Loop P|| with the No-Baryon linear model 
(from equation 1 12p . Again, since we are primarily inter- 
ested in the scale dependence of the bias, we re-normalize 
all curves to be unity at A: ~ 0.01/iMpc^^. As expected, 
the red galaxy power spectrum appears to trace the lin- 
ear theory matter fluctuations (green dash line) very well 
when k < 0.07/iMpc"^ The non-linear boost breaks this 
accordance at larger k. However, for the blue galaxies, 
the scale dependence is more complicated, having an in- 
creased pre-virialization feature and a delayed nonlinear 
boost. We also note that, for the red galaxies, the sec- 
ond and third BAOs have been almost completely sup- 
pressed, whereas only the third peak has been removed 
for the blue galaxies. The peaks and troughs, however. 
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appear to be in the right places. 

So far, we have neglected the contribution from the 
constant power 1-Halo term. In the bottom section of the 
panel we now take this into account and show Pf^ + P|| 
ratioed to the No-Baryon model. For the red galaxies, the 
agreement between the predictions and the linear theory 
that was noticed before is now broken on much larger 
scales, k ~ 0.04/iMpc"^ The trough of the first BAO 
has been shifted slightly to smaller k and the non-linear 
boost occurs at a larger scale. Because the 1-Halo term 
is about 5 times smaller for the blue galaxies, the mod- 
ifications are not as severe. The addition of this term 
offsets the suppression of power caused by the negative 
bfb2 term and the blue galaxies now appear to trace 
the linear theory on scales k < 0.07ft.Mpc^^ quite well. 
At larger k the linear spectrum is a poor match to the 
predictions. We also note that the BAOs are further sup- 
pressed and the second trough has been shifted to lower 
frequencies. 

Because PiLoop does not provide a very accurate model 
for the true non-linear power spectrum, we have studied 
the effect of exchanging PiLoop for the halof it ,54| power 
spectrum. The results are shown in the right hand panel 
of Figure [9l Although the predictions are qualitatively 
very similar, halof it predicts enhanced pre-virialization 
and smaller non-linear boosts. In the middle panel, where 
the 1-Halo term is not included, we see that the red and 
blue galaxies do not match the linear theory as well on 
large scales. In particular the blue galaxy power is sup- 
pressed on all scales except the largest. We also see that 
the BAOs have been better preserved. Although there 
are slight shifts in the positions of the second trough and 
third peak. However, the bottom panel shows that once 
the 1-Halo term has been included, the red galaxy predic- 
tions are almost as before. The blue galaxies still show a 
reasonably strong pre-virialization feature, but, because 
of the weak nonlinear boost, they match the linear theory 
rather well over nearly all the scales considered. 

We note that this modification is not entirely self- 
consistent and as such is not meant to be blindly trusted, 
since the non-linear bias terms are still derived from the 
1-Loop Halo-PT. However, we use this operation to high- 
light that a more advanced understanding of the non- 
linear power spectrum does change the results quantita- 
tively. This implies that a more advanced model of the 
scale dependence of the bias will also further modify and 
improve the predictions. 



IX. SCALE DEPENDENT BIAS AT 2 > 

We have shown that the BAO harmonic series in the 
power spectrum at {z = 0) is affected by nonlinear effects 
from bias and gravitational evolution. Although naively 
one might expect that at higher redshift nonlinear effects 
are less important, this is not necessarily so. First, one 
must pick a criterion for how to compare things at dif- 
ferent redshifts. A natural choice is to use objects of the 



same number density. In this case, as it is well-known, at 
higher redshifts objects of the same number density are 
more biased, leading to stronger nonlinear bias effects, 
even though the dark matter has less nonlinear evolu- 
tion. Therefore, overall it is not clear a priori whether 
the situation improves or not. Figure [TOl shows the halo 
power spectra at z = 1 and z = 2 measured from our 8 
LR simulations, where the halo samples were harvested 
so that they would have the same fixed comoving number 
density as the Bin 1 sample at z = (see from Table HI . 
The power spectra analysis was identical to that as de- 
scribed in Section IIIII This figure clearly demonstrates 
that the non-linear bias effects that are present at z = 
(Fig. [1]), remain present in the high redshift halo samples. 
In light of this, we anticipate that low mass halo samples 
at higher redshift, constructed so that M < M*{z > 0), 
will likewise show enhanced pre-virialization (Af*(z) is 
defined to be the halo mass at which a{M) = 1). We 
reserve further details of this issue for future work. 



X. CONCLUSIONS 

In this paper we have explored in detail, through both 
numerical and analytic means, the scale dependence of 
the nonlinear dark matter, halo and galaxy power spectra 
on very large scales k < 0.15/iMpc~^. For our numerical 
work we used and ensemble of 20 simulations in boxes 
of side 512/i~^Mpc and 8 simulations in boxes of side 
~ Ih^^Gpc. Each simulation contained more than 134 
million particles. We have found that: 

• The non-linear matter power spectrum is sup- 
pressed relative to the linear theory by 5 — 10% 
on scales between 0.05 < fc/(ft.Mpc~^) < 0.075 at 
the 2- a level. 

• The halo-dark matter cross power spectrum shows 
clearly that the bias of halo centers is non-linear on 
very large scales. The form of the non-linearity de- 
pends strongly on halo mass: for high mass haloes 
no pre-virialization suppression is seen; whereas 
there is an apparent 10% suppression of power 
relative to linear theory for lower mass haloes. 

• To make robust statements concerning their large 
scale clustering of haloes, it is essential to char- 
acterize the shot noise correction to high preci- 
sion. For high mass haloes this correction is sub- 
Poissonian, so the simple and widely used 1/n 
model must be inappropriate. Halo exclusion ef- 
fects lead to a plausible explanation for this phe- 
nomenon, which we used to motivate an alternate 
correction. However further work is required to es- 
tablish this robustly. In addition the true answer 
may need to take into account the way in which 
haloes are identified in simulations. 

• The large-scale bias of P^^ is not expected to be 
the same as that of P^^ due to nonlinear deter- 
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FIG. 10: Dark matter, halo-halo and halo-dark matter cross power spectra at z = 1 (left panel) and z — 2 (right panel) as a 
function of wavenumber. Symbol styles are as in Fig [1] At both epochs we computed the halo spectra using the same fixed 
number density of haloes as those in Bin 1 sample from {z — 0): n = 3.43 x 10~^ [/i/Mpc]^. Non-linear mass and bias corrections 
remain present at the level of 5-10% in the halo spectra. 



ministic bias; this complicates studies of stochas- 
tic bias. The difficulty of performing the halo- halo 
shot noise correction means we are unable to make 
a strong statement about the non-equivalence of b^^ 
and b^^. 

• As wavenumber increases, low mass haloes are in- 
creasingly anti-biased and high mass haloes are in- 
creasingly positively biased. Therefore, the non- 
linear scale dependence of halo bias is not simply 
due to the nonlinear evolution of the matter fluc- 
tuations. 

• Baryon acoustic oscillation features in the power 
spectrum are erased on progressively larger scales 
as halo mass is increased. In addition, small shifts 
in the positions of the higher order peaks and 
troughs occur which depend on halo mass. 

In the second half of this paper we developed a 'phys- 
ical model' to explain and reproduce these results. The 
model was constructed within the framework of the halo 
model and we focused our attention on the clustering of 
the halo centers. The halo-halo clustering term was care- 
fully propagated into the non-linear regime using T-loop' 
perturbation theory and a non-linear halo bias model. 



Our model can be summarized as follows: The density 
field of haloes was assumed to be a function of the local 
dark matter density field. Under the condition of small 
fluctuations, it was then expanded as a Taylor series in 
the dark matter density with non-linear bias coefficients 
bi{M) [l^. The density field was then evolved under 3rd 
order Eulerian perturbation theory and this provided the 
3rd order Eulerian perturbed halo density field. We then 
used the model to derive the halo-halo and halo-dark 
matter cross power spectra up to the 1-Loop level. This 
lead to the following conclusions: 

• For bi non- vanishing, the effective bias on very, very 
large scales for the halo spectra is not simply bi , but 
also depends on &2, ^3 and the variance of fluctua- 
tions on scale R. The halo center power spectrum 
contains a term that corresponds to constant power 
on very large scales. This implies that as fc — > 0, 
halo bias should diverge as [-PLin]"^^^- The halo- 
dark matter cross power spectrum does not exhibit 
this behavior. The predicted bias from this statistic 
approaches a constant value on large scales. 

• When evaluated for a realistic cosmological model, 
with non-linear bias parameters taken from the 
Sheth-Tormen mass function, the theory is in broad 
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agreement with numerical simulations for a wide 
range of halo masses. 

• The non-linear evolution of the BAOs was also ex- 
amined. The model shows that non-linear bias and 
non-linear mode-mode coupling increasingly damp 
the BAOs as halo mass is increased. In addi- 
tion, the positions of the peaks and troughs can 
be shifted by small amounts which depend on halo 
mass. 

Using the ensemble of simulations we constructed scat- 
ter plots of the halo versus dark matter over densi- 
ties, contained in top-hat spheres of size R (see Ap- 
pendix|B]). From these, it was shown that for filter scales 
R < 60/i"^Mpc the bias was indeed non-linear, and that, 
whilst the scatter increased as R decreased, the mean 
of the relationship did not change until R < 20/i~^Mpc. 
We also examined whether the nonlinear bias parameters 
derived from the Sheth & Tormen mass function [3, [l^ 
provided a reasonable match to the empirical halo bias. 

Using semi-empirical bias parameters as inputs for the 
analytic model it was shown that the model well repro- 
duced the scale dependence of the halo-dark matter cross 
power spectra in the simulations, and for all bins in halo 
mass. However, it was only qualitatively able to repro- 
duce the scale dependence of the non- linear halo bias. 

The 1-Loop halo center power spectrum was then in- 
serted into the halo model framework and defined the 
'1-Loop Halo Model'. This was used to predict the scale 
dependence of blue and red galaxy power spectra. Plausi- 
ble models for the Blue and Red Galaxy HODs were used 
and the results showed complicated scale dependence. 

Significant work still remains to be performed for this 
analytic approach to be sharpened into a tool for pre- 
cision cosmology. Some possible improvements are: ex- 
changing the 1-Loop matter power spectrum for an ac- 
curate analytic fitting formula i.e. after the fashion of 
halofit, but designed purely for large scales; an ap- 



plication of the new re-normalized perturbation theory 
techniques [4l|, [o^] coupled with the non-linear bias 
model should certainly produce better results. 

The analytic and numerical results that we have devel- 
oped here are concerned purely with the clustering in real 
space. In a subsequent paper we shall extend our analysis 
to explore the more realistic situation of the scale depen- 
dence of dark matter, halo and galaxy power spectra in 
redshift space. 



When this paper was nearing completion, a preprint 
appeared (Qlj with similar calculations of how nonlinear 
bias changes the power spectrum. 



Acknowledgements 

RES would like to thank Peter Schneider for a very 
thought provoking discussion that has led to the current 
paper. We thank Jacek Guzik, Gary Bernstein, Bhu- 
vnesh Jain and Laura Marian for many useful discussions 
during this work. We would also like to thank Pat Mc- 
Donald, Daniel Eisenstein, Eiichiro Komatsu and Martin 
White for useful comments on the draft. RES thanks 
Joerg Colberg and the Virgo Consortium for providing 
access to the Hubble Volume and the VLS simulations. 
RES thanks CCPP and NYU for their kind hospitality 
during part of this research. RS would like to thank 
M. Manera for useful discussions on halo exclusion. RKS 
would like to thank Martin White for many illuminat- 
ing discussions and the Aspen Center for Physics. We 
thank M. Crocce and S. Pueblas for help regarding the 
numerical simulations used here. RES and RKS both 
acknowledge support from the National Science Founda- 
tion under Grant No. 0520647. RS is partially supported 
by NSF AST-0607747 and NASA NNG06GH21G. 



J 



1. Poisson model 

It is necessary to correct the measured dark matter power spectra for shot-noise errors. These arise through 
approximating the continuous CDM fluid by a point process. If the discretization of the density field obeys a Poisson 
process, i.e., dark matter particles are placed with probability oc p[l + S{x)]SV, then the 'discreteness' correction is 
|: 

^'truo(fc) -Pobs(^) - -Pshot ; P'shot = l/n . (73) 

Since the number of particles in our simulations is large, 512"^, this correction, on the scales of interest, is insignificant, 
e.g. Pshot/-P(0.1/iMpc~^) ~ 10"'^. However, one must also correct P^^^ for discreteness. If the dark matter haloes are 
also regarded as a Poisson sampling of the smoothed halo density field, then the correction will be the same but using 
the appropriate number density fih- In Figure [Til we show the effect of this standard correction on our halo power 
spectra and we plot the ensemble of the shot-noise corrected halo power. The negative power values that result at 
high k demonstrate that this model must not be exactly correct and therefore for accurate measurements it must be 
modified in some way. We now discuss a possible explanation for this and propose a new shot correction. 
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FIG. 11: The effect of Poisson Shot noise correction on the halo-halo power spectra for the four bins in halo mass. 
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APPENDIX A: HALO DISCRETENESS CORRECTIONS 



1. Halo exclusion effects 



As we have seen Poisson sampling must not be exact for haloes, particularly those of large mass. A possible reason 
is that by using a friends of friends algorithm, one is automatically imposing that haloes are never separated by 
distances smaller than about the sum of their radii, else they would have been linked as a bigger halo. As a result of 
this exclusion effect, the two-point correlation function of haloes drops dramatically from a value ^ 1 at r ~ re to 

= — 1 for r < re, where Ve is the exclusion radius. This sharp drop of has an impact on the large-scale power 
spectrum, as we now show. 

Absent exclusion, the power spectrum would have been 



However, exclusion effects mean that 



kr 



— Wdr-i;TH(fc), 
kr 



(Al) 



(A2) 
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where VTnik) = (47r/3)rg WTH(fc?'e) with Wth the Fourier transform of a top-hat window in real space. Hence, the 
difference in power due to exclusion is 

Apf^ik) = r ^^4^ (1 + C") 47rr2 dr. (A3) 
Jo kr 

To estimate this, we must model <^'* at scales smaller than the exclusion imposed by the friends of friends definition 
of a halo. We do so simply by approximating ^'^ by a power-law ^'^ = £,e{r/re)~^ obtained from fitting to the r > 
measurements, where is the exclusion scale. In this approximation 

AP''{k) = VTii{k) + v{k) ^e, (A4) 

where, in the large-scale approximation kr^ <C 1, 

vik):.^--^(kr.f. (A5) 
3-7 3(5-7) 

Setting 7 = gives wth in the same limit. Thus, we see that exclusion makes the halo power spectrum smaller by 
the amount given by Eq. (jA4p . In contrast, Poisson shot noise makes the power larger by 1/nh. If Poisson noise is 
subtracted from halo power spectra in simulations, then the result becomes negative at high k. So it may be that 
adding back the power lost to exclusion will make the power positive again. This is our procedure. Note that the 
noise in Eq. (jA4j is not white. 

As a final note we reemphasize that all of these troubling issues may, to a certain extent, be neatly side stepped, if 
we measure the halo-dark matter cross power spectrum. The clear advantage of this approach is that the 'shot-noise' 
is dramatically reduced (there are many more particles than haloes) and exclusion no longer plays a role. 



APPENDIX B: MEASURING THE NON-LINEAR BIAS PARAMETERS 



The accuracy of 1-Loop Halo PT depends on the accuracy of the halo bias parameters bi, 62 and 63. Figure |4] shows 
these parameters derived from the conditional Sheth & Tormen mass function [l5|. There are a number of ways to 
test the accuracy of the bias parameters. The most direct is to smooth the halo and mass fields with a filter of scale 
R, and to then plot S^{x\R) versus S{x\R) in the smoothed fields (e.g. [zEUll)- From the resulting scatter plot, one 
then directly fits for the parameters bi up to some order, subject to the constraint that (J*^) = (S) = 0. 

Figure [12] shows such scatter plots for the four bins in halo mass described earlier (Table HI. In all cases, unlike 
elsewhere in the paper, the fields were filtered with a 'top-hat' window function. Solid blue, red, green and magenta 
lines show results for four filter scales Rtr = {20, 30, 40, 50}/i~"'^Mpc. The upper and lower dotted lines show the 
1 — (T errors about the mean relation (solid lines). (Smoothing with a Gaussian filter instead does not significantly 
change the mean, but does slightly change the scatter around the mean relation.) The figure shows that, above a 
certain scale, (^'^l^) defines a universal curve which is independent of R; it is this curve which the bias coefficients 
are supposed to describe. The figure also shows that the scatter around this mean curve decreases as R increases. 
The assumption that bias is deterministic is equivalent to assuming that this scatter, this stochasticity, is negligible. 
Clearly, this is a reasonable approximation only for large smoothing scales R, and R must be bigger for the rarer, 
more massive haloes before the bias can be called deterministic. 

The main problem with this approach is that, when R is large, then the distribution of S becomes sharply peaked 
about its mean value (zero!), so estimating the higher order bias coefficients becomes difficult — only the linear bias 
term can be reliably measured. In effect, determination of the bias coefficients corresponds to fitting a polynomial 
to these mean curves, and the best-fitting coefficients will be correlated. The results of this exercise are tabulated 
in HTl which also shows the values of bi predicted from the Sheth- Tormen mass function. Some of the discrepancy is 
a consequence of the fact that the Sheth- Tormen mass function under-predicts the abundance of high mass haloes 
in the HR and LR simulations by up to 20% (see 52] for an explanation as to why this happens). But note that, 
because of the correlation between the fitted coefficients, the discrepancy with the theory values is difficult to assess. 
The solid thick (cyan) line in the figure shows the mean bias relation derived from the Sheth- Tormen mass function. 

In light of these issues we use the following prescription when generating the theory curves: 61 coefficients are to be 
measured directly from the data as described in Section [I IIBI and the parameters 62 and 63 are to be derived through 
computing the following integrals over the relevant bin width: 

= -(M ^^^ / dMn{M)h{M) , (Bl) 
n{Mi,M2) Jm, 
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FIG. 12: Scatter plot of halo number-over-density vs the over density of dark matter for four bins in halo mass. The fluctuations 
were measured in top-hat spheres of radii Rth ~ {20, 30, 40, 50}/i~^Mpc. The means and 1-sigma errors measured from the 
ensemble are denoted by the solid blue, dash red, dot-dash green and dotted magenta hnes, respectively. The solid thick cyan 
line shows the predictions from the Sheth & Tormen model. The thick triple dot-dash lines show the results from the cubic fit 
in S as described in the text. 



where n(M) and bi(M) are the Sheth & Tormen mass function and bias parameters. We shaU refer to this approach 
as the semi-empirical method. It shall be left to our future work to to provide a more self-consistent solution to this 
problem. 



APPENDIX C: EULERIAN PT 



1. Real space kernels 



The growth of density inhomogeneities in an expanding universe may be explored in the single stream approximation 
using the Eulerian fluid equations. In Fourier space the equations governing the evolution of the fluctuations to the 
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density (5(k, t)) and divergence of the peculiar velocity field (6'(k, t) = V • v) are written |10 |: 
dS{k, t) 



dt 
d0(k,t) 
dt 



(27r)3 



[(5^(k)]2«(qi,q2)e(qi,i)<5(q2,t) 



H{a)0{k,t) + ^n„M\t)6{k,t) = - J [<5^(k)]2/3(qi,q2)(?(qi,<Mq2,i) 



where 



aij = a(qi,qj) 



k • qi 



; Aj=/3(qi,qj 



(27r)3 

fc^(q» • qj 



(Cl) 



(C2) 



(C3) 



As shown in [lO|, for the case of the Einstein-de Sitter model the density and velocity divergence fields may be 
expanded as a perturbation series and the solutions at each order may then be written down in terms of the density 
and velocity perturbations from all lower orders: e.g. 



<5(k) = ^i?J'(a)<5„(k) 



<S„(k) 



n=l 



{2n) 



3n — 3 



[<5^(k)]„F(/)(ki,. 



0(k) =i/(a)/(a)f;n,"(a)^?„(k) ; &„(k) = / ^^-'fl'^'^lf'^^ [<5^(k)] „ G,(^) (k^, k„ 

n=l ^ ' 



(C4) 
(C5) 



where H(a) is the dimensionless Hubble parameter, Di{a) = a is the linear growth factor and /(a) = d log Z^/d log a « 
i7m^ gives the velocity growth suppression factor. The functions Fn^ and Gn^ are the PT kernels for 5 and 9 
symmetrized in all of their arguments, respectively. The notation [(5'^(k)]^^ = (5^(k — qi — • • • — qn) was adopted. 
The first three symmetrized density kernels are: 



FW(q) = l; 
i^i^^(qi,q2) = 
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[ai_2 + ^24] + 3/3i,2 ; 



i^i'''(qi,q2,q3) = 

2 r 

"27 



7' 

i^i':^a3,i2 + J^i'i«i,23 + i^3aa2,3i + G^^lai2,z + G^'^a23,i + 4!la3i,2 



^0/^12, 3 + G'2'^3/323,1 + G'3'^1/33 



The first three symmetrized velocity divergence PT kernels are: 
Gi^^(q)-1; 



3 4 
G2'^(qi,q2) = [ai,2 + a2,i\ + y/3i,2 



4'^(qi, q2, q3) - ^ -^^2^3,12 + J^2>1.23 + i^3"l«2,31 + G'(;^ai2 



1 
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2 

9 



•-^2, 3*^23,1 + (-^3, 1031, 2 



^1*^2/312, 3 + G2\P2i,l + 03*^1/331,2 



where we have adopted the compact notation 



Fi^)(q,,, . . . ,q,J ; G^_,^ ^ G(f^(qn, • ■ • , J 



and 



Hi...ln,]l---3r, 



a(q«i H ^q*„,qji h ^qj„.) 



in-.j = /3(q,:i H ^qj„,qj) 



(C6) 
(C7) 

(C8) 

(C9) 
(CIO) 

(Cll) 

(C12) 
(C13) 



Exact analytic solutions for arbitrary cosmological models have not yet been found. However, as was shown in [93| . 
under the assumption that Dn(a) oc [Di(a)]" and ~ ^ solutions are identical to the Einstein-de Sitter 

solutions, but the growth factors are changed to those for the particular cosmological model in question. All other 
changes are small corrections [loj . 
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2. Evaluation of 1-Loop expressions 



Some care is required in the numerical evaluation of the Halo-PT expressions because fine cancellations can occur 
between negative and positive terms. The approach that we adopt throughout this paper can be demonstrated through 
the following example. Consider the 1-Loop power spectrum for CDM given by equation The PT corrections 

Pi3 and P22 maybe analytically developed up to the following points: 



P22(fc) 



252(27r)2 

(27r)3 
°° dq 



r dxx^Pnixk) 1-42x2 + 100 - 1^ + + -(1 - x^)^7x^ + 2) log [^^1 1 

Jo { x^ x'^ X Lf-i|JJ 



(C14) 



(27r)2 



Pii(g)Pii(|k-q|)[P2(q,k-q)]^ 



dqPii{q) / d^Pii(fc-0(a;,^)) 



5 1 ^ — X 
7^ 2V(x,/i) 



■0(a;,M) 



ipix,fi) 



2 


/i — a; 




.+7 


_^p{x,fi)_ 





,(C15) 



where x — q/k and where 'ip'^{x, ^) ~ 1 + x'^ ~ 2x^. If the P13 integral is truncated on large and small scales to inhibit 
infra-red and ultra-violet divergences, as may occur for some power spectra, then identical constraints must also be 
placed on the P22{k). Explicitly, if we adopt 



P(q) = for 



k < fcfun 
k > fccut 



then the angular integral for P22 must necessarily have the new limits 



IJL2 — mm 



/ii — max 



(C16) 



Lastly, the particular values for the limits fcfun and fccut were selected so that variance integrals, a^{R), would be 
convergent within the finite range. 
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1 I |2 

1=1 

where M are the number of Fourier modes in a spherical 



shell in fc-space of thickness Afc, and * denotes complex 

conjugation. 

[96] The procedure was also performed separately for the no 
baryon model spectrum, this was done in order for ratios 
to be taken. 



